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In this paper we make use of the concept of conserved geometric phase and of group velocity,in con- 
junction with the representation theorv[44ll4{T|. in order to derive some relevant physical quantities 
for the description of the dielectric and magnetic response of crystalline materials. As an application 
of the model, we derive the expression of the macroscopic dipole moment per unit volume, and the 
expression of the current induced by a uniform static external electromagnetic field. 



INTRODUCTION 

In 1983, in a milestone paper, M.V. Berry|l| showed 
that a quantal system in an eigenstate, slowly trans- 
ported round a circuit C, by varying parameters f3 in its 
Hamiltonian Hp, will aquire a geometrical phase factor 
$ in addition to the dynamical phase factor. 

Moreover, in 1986 Y. Aharonov and J. AnandanQ 
showed that it is possible to define a new geometri- 
cal phase factor for any cyclic evolution of a quantum 
system, independently of the adiabatic approximation Q 
and expressible as a gauge invariant quantity, including 
Berry's case as a particular case. 

In 1988, J.Zak 3] showed also that it is possible to de- 
fine a Berry's phase for the dynamics of electrons in infi- 
nite periodic systems, suggesting how to label bands of a 
crystal by making use of the conserved geometric phase. 

In early nineties R.Resta and D.Vanderbilt 0, 
showed, how to relate the conserved geometric phase to 
the macroscopic dipole moment per unit volume of a crys- 
talline insulator. The model have been called by the au- 
thors the modern theory of polarization. In the following 
section, we shall review the concept of polarization in 
crystals, within the modern point of view 
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POLARIZATION AND GEOMETRIC PHASE 

The problem of the dielectric and magnetic response of 
matter, in condensed matter physics, have been inten- 
sively studied in the last decades (and references 
therein). A particularly interesting example is the treat- 
ment of the electronic structure of a crystalline material 
in an external uniform electrostatic field. The modern 
theory of polarization in crystals^, Q is one of the most 
important recent developments, in fact, the macroscopic 
polarization is a basic quantity used to describe, classi- 
cally, dielectric media. 

However a crystal is only a particular example of a 
quantal body and the classical concept of polarization 
extends also to non crystalline (disordered) systems, and 
it should not depend on the particular model system em- 
ployed to perform explicit calculations. 



As shown by Restafly it is possible to define variations 
AP, of the dipole per unit volume of a crystal, due to a 
source, and given by 



AP = P (1) - P (0) 



f 1 dP 



(1) 



Here A is a dimensionless parameter depending on time; 
A = and A = 1 represent the initial and the final states 
of the system, along a transformation that changes the 
periodic potential. 

An explicit representation of p( A ) in terms of the am- 
plitudes u„k of Bloch-statesj23,0,|^, has been derived 
by King-Smith and Vanderbilt |fjj: 



(27T) 



3 A^i 



dk<u^V k \u^> (2) 



(A) 



BZ 



where e > is the absolute value of the charge of the 
electron, BZ is the Brillouin zone of the crystal and n is 
a discrete index, running over occupied bands. 
Variations of this dipole per unit volume have been in- 
terpreted as the integrated electronic polarisation cur- 
rent 0|. According to the modern point of view, 
only variations of polarization are physical observables, 
completely independent of the periodic charge distribu- 
tion. The charge density is represented by the squared 
modulus of the wave functions 4] . Variations of polariza- 
tion, instead, have been related to a Berry's phase^, Q 
of the system. We shall try to clarify this statement in 
what follows. 



GEOMETRIC PHASE AND V MATRIX 

Here it is shown how to develop a calculation scheme 
in order to model the response of a body, in its 
steady states [3lj|. to a static uniform electromagnetic 
field. The model is based on the concepts of group ve- 
locity and of conserved geometric phase. 

At this point we point out that we may assume defini- 
tion JQ| of being valid even when a uniform static external 
field, is present in the system, as it is independent, i.e. 



2 



external, to the charge distribution of the quantal body 
(QB). In the following we shall always consider the states 
of the body of being implicitly dependent on external 
fields. 

It is already known, in condensed matter physics, that 
the macroscopic electronic dipole moment per unit vol- 
ume, of an insulator is related to a conserved geometric 
phased 

Here, we show how to obtain a conserved ge- 
ometric phase making use of a matrix that we shall call 
V matrix. V is obtainable from the expectation values 
of the group velocity operator [si). 

We start our discussion presenting in this section part 
of the work of ZakjjJQ], in order to introduce the rea- 
soning which led us to dervive results reported in the 
following section. Firstly, we derive the conserved geo- 
metric phase in absence of external electromagnetic fields 
in the system, being the system an infinite crystal. Let 
us start writing the Schrodinger equation in the k — q 
representation 0,0 as follows 



and because of the normalization condition of the com- 
plex wave functions u n k follows that 



d* 



d nr , 



(8) 



also implying that the diagonal fields are real. From 
eq.© follows 



m 7^ n 



(9) 



h 



Let us first note that the trace of the velocity operator, 
evaluated with respect to the occupied bands, is 



Tr(v mn ) = ^v nn = ^2 *rVk£n = v(k) 



(10) 



[— (-iRV, + Kk) 2 + V(q)]u„, k (q) = H 0M u n ^(q) 



2m 



£ n iVk(q) 



(3) 



In analogy to Q, , if we now derive eq. (JSJ with re- 
spect to k we obtain 



— (-iSV q + 7jk)tt n ,k(q) + ffo,kV k w„,k(q) = V k e„M„,k(q) 

+ e«V k u„ ik (q) 
(4) 

Multiplying eq. J3J on the left by , dividing by h and 
integrating with respect to q in the super-cell we obtain 



that is an invariant quantity because of its geometric 
nature. It represents the electronic current density of 
the system, in k space, once multiplied by the electronic 
charge e, 



j(k) = eTr(v mn ) = ev(k) 



(11) 



By relation (|45|l in the appendix, and taking the curl 
of vectors d nn it si possible to define 



V k x d nn = i J V k < k (q) x V k u„, k (q)dq (12) 

as already shown by Zak^ij. Also, by equation JJJ 
follows that 



Vmn =< Um,k|v|w n , k >= -lLO mtl d mn + -V k e„<5 mn (5) 

n 

where we write the transition frequencies as Lu mn = 
; '"7 £ " and 



iV k u„ ik (q) = y^d mn u m , k (q) 

m 

and combined with ea. l|12(l gives 



(13) 



{-ihVq + Kk) 



rn 



(6) 



is the group velocity operator as shown in ref.|3l|. 
Coefficients v mn are the matrix elements of the V ma- 
trix. 

In the following, we shall express < u m ,k|v| , u nik > as 
also < m|v|n > in order to simplify the notation. 
Matrix elements d„, n are defined as 



: < Um,k|«Vk|Wn,k > 



(7) 



V k x d nn = i ^2 d mn x d r 

Making use of eq.© we obtain 



V k x d nn = i d mn x d nm = i 2 



(14) 



(15) 



in agreement with Berry and Zak's derivations 0, • 
We can evaluate then the flux of the field given by 

ea. (|14fl . through the closed surface S of the Brillouin zone 

in k space, containing a closed curve C. 
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We may either write 



I rf S .Im^[ Vmn ^ V " m ]= f dS-V k xd nn 

= <j> d nn • d\ = <j) n {C) 



(16) 



or in a different way, 



J dS ■ Im ^2 d mn x d nm = J dS • V k x d nn 

= j> d nn • d\ = (f> n {C) 

(17) 

where Im stands for imaginary. After a summation 
over n, we find 



^ n (C) = <I>(C) 



(18) 



meaning of macroscopic dipole moment per unit volume 
of the system, in its steady states. 

Berry's phase instead is given by eq.J2U, so that we 
may think of a classification of materials, as a function 
of values of $(C), directly related to definition i|20|). 

Moreover, we may calculate the conserved phase &(C), 
not only by the knowledge of matrix elements of V, as 
in eq. I|16|) , but also by the knowledge of matrix elements 
d mn by eq. fTTjl . A discussion about degeneracies of the 
energy spectrum can be found in references 0, 0, E3- E~l 

El- 

For what stated in 0, and by results given in [3l| 
the reasoning may be directly extended to a many-body 
theory. 

In that context the expression will be generalized 
by eq.(18) of HU- 



v m „ =< m|^V k iJ|?i > 



(22) 



being H the k dependent Hamiltonian of the problem. 

In the following, we extend the problem, to the case 
of an external uniform magnetostatic field present in the 
system. 



being $(C) is the conserved geometric phase. In view 
of eq. (|15f) we note that V k x d nn is the trace of the Berry's 
curvature^, while instead d nn is Berry's connection [J. 
Let us calculate the averaged trace of fields d nn , inte- 
grated over the Brillouin zone in reciprocal space, or in 
other words 



V MATRIX IN MAGNETIC FIELDS 

In this section we generalize the problem to the case 
of an external uniform magnetostatic field B° present in 
the system. The Hamiltonian operator is given by 



P = ed 



OCC „ 

k^L dk 



< u„ k |iV k |u„ k > 



dk 



(19) 



= e / - , d 

IBZ 



(2tt) 3 



where e is the electronic charge. The circuit integral 
of the field P defined as 



^.[Ihw^^Vx^ (23) 

+ F(q)] 

The corresponding expression of the group velocity op- 
erator is then, 



V (-B°) = I Vk tf k 

n 



(24) 



so that matrix elements of the velocity operator are 
expressible as follows 



(20) 



IS 



,(B U 



=< m v v ; \n > 



1. 



-V k e n <5 n 



(25) 



j P • dl = e$(C) 



(21) 



Ea. H19(l is coincident with results found in ref.(0,0|) 
(see also 0, 113, 13 f° r ^ s use m numerical cal- 
culations). Because of that we can associate to P the 



— B° x d r 



In analogy with (fTTJl , the current density is 



j(k) = eTr(v mn ) = ev(k) + -oj cycl b x P (26) 
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being b a unit vector parallel to the direction of the 

e\B°\ 

magnetic field and uj cyc i — — the cyclotron frequency 
of an electron. Moreover, by equation (|26[1 and l|48|l we 
can calculate the divergence of the current as follows 



V • j(k) = eV • v - -u> cyc ib ■ rotP 



(27) 



In the following section we shall derive another vecto- 
rial field, that is directly related to the curl of the dipolar 
field H20f> . For the sake of clarity we report in the text 
only main results, while details about our calculations 
are given in the appendix. 



we find 



= —B° x rot d„ 

mc 



,(B°)1 - 



[d n 



X VJ, 



(34) 



directly derivable by d mn and v mn matrix elements, 
it is clear that 1^ = and we can define 



(35) 



J s FIELD 



Let us define 



(B°) (B°) (B°) (B°) 
Vmn X V nm -)- V nm X V mn 



2u)r 



(28) 



where S stands for symmetric combination. After a 
summation over ra/nwe obtain 



occ (B°) (B°) (B°) (B°) 
^""^ jV mn X v nm + v nm X v mn , 

e 



(29) 



B° x rot d 



mc 



nn 



being rot d nn = V k x d nn . 

By eq. <|19|) is valid the following relation 



-b x rot P = J s 



where 



(30) 



(31) 



Definition l|3U|l also implies that 



dVV ■ J s = b • / dS x rot P (32) 
Js 

We associate the meaning of a current to the field J s , 
whose source is related to the rotational properties of 
P. An alternative derivation of the above results can be 
obtained by the definition of another vectorial field. By 
the following relation 



[Vmn" 1 X d nm] = -^mn(d mn X d nm ) 



mc 



-(B° x d mn ) x d r 



(33) 



APPEARANCE OF THE GEOMETRIC PHASE 

Taking the vectorial product of B° and the field we 
obtain 



B» x 



e\B 



0|2 



e\B 



mc 

0|2 



mc 



-b x (b x rotd nn ) 
-(b • rot d nn )b - 



mc 



-rot d r 



(36) 



Let us exclude the cases rot P = 0, and rot P parallel 
to b. 

Dividing ea. (|3(jH by summing over n and calcu- 

lating the flux of over the surface S of the Brillouin 
zone, containing a closed curve C we obtain 



-[/ dS-(bxv s )] = ^ / 
Js c Js 



dS ■ rot d 



^cycl 

c J„ 

^cycl 
C 

^cycl 
C J,-, 



[ (b ■ rot d)b • dS 

Js 

/ (b • rot d)b • dS 

Js 



(37) 



Multiplying eq. IpTfjl by the electronic charge and bear- 
ing in mind ea. H21j) we obtain 



b • / dVrot J s = -{[ dS ■ (b x J s )] 
Jv Js 

-e$(C) 

(b • rot P)b • dS 



^cycl 
C 

^cycl 



(38) 
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LOW POWERS OF ± 

c 

Bearing in mind definition l|25|l . we can write at first 
order of powers of - 



Defining 



^™^d mn x d r 



e (39) 

+ ^mn -B X (d mn X d nm ) 

mc 



(B°) (B°) (B°) (B°) 
r V mn X v nm v nm X v mn 



(40) 



dividing by w^ n , and summing over m ^ n we obtain 



= rot d nn 



occ (B°) (B°) _ (B°) (B°) 
Vmn X V nm V nm X V mn 



9/ >2 
AUJ mn 



(41) 



Evaluating the flux over the surface S of the Brillouin 
zone we find 



/ dS ■ [V£] - <M<7) 



(42) 



If we instead consider also terms containing second 
powers of -, we instead obtain 



dS ■ [V*] 



c 2 ^ JS 



rfS-[-f^(b.(d mn xd n 



0)b] 

(43) 



Equations and (|38[1 . in conjunction with ea. i|21|) 
and ea. l|3U|) . clearly show that if the vectorial field d is 
conservative (irrotational, rot d = 0), then the system 
would acquire a vanishing Berry's phase, and the vecto- 
rial field J s would become vanishing by ea. H3U|) . More- 
over, if d is conservative, the divergence of the current 
j(k) will not be affected by the latter. In the case of d 
not being conservative the flux of j(k) over the surface of 
the Brillouin zone will be instead proportional to the pro- 
jection, over the direction of the applied magnetostatic 
field, of the circulation of d (or equivalently P) over the 
surface of the Brillouin zone. 



POWER IN A UNIFORM STATIC 
ELECTROMAGNETIC FIELD 

The group velocity is not explicitly dependent on a 
uniform electrostatic field Eq eventually present in the 



system, so that all the results obtained above are equally 
right even if Eq is present in the system. Bearing in mind 
ea. (|26[) and ea. H19l) . we can express the work per unit 
time and per unit volume done by the electromagnetic 
field (E , B°) on the body as 



W = E ■ J 

= eE° 



E° 

V 
dVv(k 



dVj(k) 

47re 2 . 



dVd 



(44) 



47re" 



S • d 



being S = ^(E x B°) the Poynting vector [H Eg 
associated to the electromagnetic field and V the volume 
of the Brillouin zone in /c-space. 



CONCLUSIONS 



The reported non-relativistic calculation scheme, may 
be useful for practical calculations of macroscopic electro- 
magnetic properties of a quantal body interacting with a 
static external uniform electromagnetic field. It has been 
shown how to calculate the power density, the macro- 
scopic polarization P, and the electronic current induced 
by the presence of a uniform magnetostatic field B°, mak- 
ing use of the concept of group velocity. It has been 
shown the reasoning in the case of the one-body approxi- 
mation, and by construction it would apply to every first- 
principles calculation scheme 0, IS H> Ell > based on peri- 
odic boundary conditions within a DFT scheme[3(J. The 
above discussion can be extended to a 'many-body' the- 
ory because of what stated in[||. Ea. (|16[) may help us 
for a theoretical classification of crystalline materials via 
the values of $(C), that are directly dependent on the 
dielectric properties of the material under consideration 
via eq.J2U- We may also think that dielectric and mag- 
netic properties of the crystalline material determine its 
geometrical properties too. In fact ea. H21() .ea. H30(l show 
that if the vectorial field d is conservative (irrotational, 
rot d = 0), then the system would acquire a vanishing 
Berry's phase and the vectorial field J s would become 
vanishing by ea. H3U|) . Nevertheless a proper theoretical 
solution of the problem requires a theoretical study of 
the ionic response as well as the electronic one, that we 
leave to a future work. The model presented above al- 
lows us to find several useful physical quantities in a clear 
and simple way, amenable of a straightforward physical 
interpretation. 
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APPENDIX 



we can recast ea.(|51jl. as follows 



Vector Identities 



(B°) (B°) _ 2 j d 
v mn x v nm — u mn u mn X Ud 



V x (/A) = /V x A + (V/) x A 



(45) 



+ —r^[{B° x d mn ) x (B° x d nm )] 



A x (B x C) = (A • C)B - (A • B)C (46) 



(AxB)x(CxD) = ((AxB)-D)C-((AxB)-C)D (47) 



(V-AxB) = B- VxA-A-VxB (48) 



iu mn — B° x (d mn x d nm ) 

mc 

" ^mn^ mn X ^nm 



:P°xd mn ).d nm )B° 



- ((B° x d mn ) • B°)d r 



-B° x (d mn x d nm ) 



+ 4rP°xd mn )-d nm )B 



(53) 



J 6 FIELD 



Let us evaluate the following quantity 



-B° x (d mn x d nm ) 



tl/ d nln X dj) 



;(B°-(d m „xd nm p° 



(B°) (B°) 2 j j 

v mn x v nm _ ^ mn u mn x U nm 

e 2 

+ -^[(B° X d mn ) X (B° X dnm)] 

(49) 

- !w mn [d mn x — (B°x d nm )] 
mc 

+ wj nm [d nm x — (B° x d mn )] 

mc 

From ea. (|49|l . making use of the vector identity 1)460. 
and the fact that w mn = -w n m we can recast the last 
two terms of ea. H49|l as follows 



i^mn[d m „ X — (B° X d nm )] 

mc 

w„ m [d nm x — (B° x d mn )] 

mc 

■ ^mn B X (d mn X d nm ) 

mc 



(50) 



Wmn B° X (d mn X d nm ) 

mc 



and define the following quantity, 



(B°) (B°) 
^ inn X V n m , , 

j > - = rot d„ 



EVr 



mn 



Y] — ^-b x (d m „ x d„ m ) 



m^n 



V oT' (b • (d mn x d nm ))b 

(54) 



m^n mn 



so that 



V (B°) x V (B°) 
mn nm 



: U3, 



2 d„ 



m"mn x d nm 



B° X (d mn X d„ m ) 

mc 

;p°xd ra „)x(B°xd nm )] 



(51) 



Making use of the vector identity (|47J) and 



(A x B) • A = 



(52) 



being b a unit vector in the direction of the magnetic 
axis and u cyc i is the electronic cyclotron frequency per 
unit magnetic field, defined as jg£y = ^ ■ 

We can write 



M B °> x v< B °) 

L mn nm 



2 j 



x d r 



+ «Wmn #° X (d mn X d„ m ) 

mc 

+ ~^(B° • (dmn X d„m))B° 



(55) 
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and 



We find 



r (B°) (B°)i _ 2 d H 
l v nm x v mn J — w nm a nm x a : 



mn 



+ «w nm — B u x (d nm x d mn ) 

mc 

+ 4^(5°-(d nm xd mn ))B° 

+ iw mn — B° x (d mn x d nm ) 

mc 

-^(B°-(d m „xd„ m ))S» 



r(B°) (B°)l 
L mn rim J 



2ia; n 



-B° x (d mn x d nm ) 



ea. l|55|l and ea. (|56[) imply then 



(B°) (B°) (B°) (B°) 
Vmn X V nm -)- V nm X V mn 

2w m n 



= « — -8° x (d mn x d nm ) 

mc 

Also, summing over m ^= n we find 



(56) 



(57) 



* I] [ V mn 0) >< dnm] - [v^ X d mn ] = — B° X rot d nn 

(63) 

In fact by making use of definition (|14p. the vector 
identity l|46|l . eq.©, and the relation w mn = — tj nm , we 
find 



[v^„ 0) x d nm ] - [v^ x d mn ] = 

= -Mmn(dmn X d nm ) + — (B° X d mn ) X d nm 

mc 

- H^nm(d„ m x d mn ) + — (B° x d nm ) x d mn ] 

mc 

- "! w mn(dmn x d nm ) H (B X d mn ) X d nm 

mc 

- [iw mn (d nm x d mn ) + — (B° x d nm ) x d mn ] 

mc 



- , mn(dmn X d nm ) 



-{B° x d mn ) x d r 



[-iWmn(dmn X d nm ) H (B° X d nm ) X d mn ] 

mc 

-iw mn (d mn x d nm ) + — (B° x d mn ) x d nm 



(B°) (B°) (B°) (B°) 
mn X V nm T" V nm X V mn . 



e 



2lj t 



= — B° x i y~] (d mn x d nm ) = — B° x rotd r 
mc *—f mc 

eq.JSSJ and ea. l(T^|l implies 



where 



Jcycl 



b x rotP = 



(58) 



(59) 



^mnfdmn X d nm ) (B X d nm ) X d mn 

mc 

: — (B° x d mn ) x d nm (B° x d nm ) x d mn 

mc mc 

: — hd„m x (B° x d mn ) + d mn x (B° x d nm )] 

mc 

^ —[-(d„m • d mn )£ + (d nm ■ B°)d mn 

mc 

(d mn • d nm )B — (d mn • B )d 

nmj 

= —[(d„m • S°)d mn - (d mn • B°)d nm ] 

mc 

: —5° x (d mn x d nm ) 

mc 



(64) 



P = ed = e ^2 dnn 



(60) 



By definition (|25|l ane eq.© follows that 



and 



s -ev s 



E X V mn° } ] - [dmn X V n *°>] 



(61) 

is a pure imaginary field. 

Also we may re- write ea. (!63|l as 



(65) 



1 A field 



Let us write the following quantity 



v mn ^ X d nm ] — *^mn(dmn X d nm ) 



mc 



-{B° x d mn ) x d r 



— B°x rot d r 

mc 



(66) 



(62) We can evaluate the following flux over the surface of 

BZ (Brillouin zone) in k space, containing a closed curve 
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C. 



n J S 

= ^ Ucyd J dS ■ (b x rot d nn ) 

= y~] uj cyc i / dW ■ (b x rot d nn ) 

= - ^ Ucycl / dvb ■ (V x rot d nn ) 

= —Ucyci / dVh ■ (V x rot d) 
Jv 

= —LUcycih ■ / dV(V x rot d) 
Jv 



(67) 



= — LU C ycib ■ / dS x rot d 
Js 

where in the third line has been used Gauss's theorem 
(divergence theorem), and u> cyc i = ^- is the cyclotron 
frequency per unit magnetic field and B = \B |b; being 
b the direction of the magnetic axis. Equations Q59JI and 
(loTjl imply also 



dW ■ J & = b • / dS x rot P (68) 
Is 
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